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Collisions between galaxy clusters provide a test of the non-gravitational forces 
acting on dark matter. Dark matter’s lack of deceleration in the ‘bullet cluster 
collision’ constrained its self-interaction cross-section < 1 25cm^/g 

(68% confidence limit) for long-ranged forces. Using the Chandra and Hubble 
Space Telescopes we have now observed 72 collisions, including both ‘major’ 
and ‘minor’ mergers. Combining these measurements statistically, we detect 
the existence of dark mass at 7.6(T significance. The position of the dark mass 
has remained closely aligned within 5.8±8.2 kpc of associated stars: implying a 
self-interaction cross-section <0.47 cm^/g (95% CL) and disfavoring 

some proposed extensions to the standard model. 

Many independent lines of evidence now suggest that most of the matter in the Universe is 
in a form outside the standard model of particle physics. A phenomenological model for cold 
dark matter ( 1 ) has proved hugely successful on cosmological scales, where its gravitational 
influence dominates the formation and growth of cosmic structure. However, there are several 
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challenges on smaller scales: the model incorrectly predicts individual galaxy clusters to have 
more centrally concentrated density profiles (2), larger amounts of substructure ( 3 , 4 ), and the 
Milky Way to have more satellites able to produce stars (5) than are observed. These incon¬ 
sistencies could be resolved through astrophysical processes ( 6 ), or if dark matter particles are 
either warm (7) or self-interact with cross-section 0.1 < < 1 cm^/g ( 8 - 10 ). Follow¬ 
ing ( 11 ), we define the momentum transfer per unit mass integrating over all scattering 

angles and assuming that individual dark matter particles are indistinguishable. 

Self-interaction within a hidden dark sector is a generic consequence of some extensions 
to the standard model. For example, models of mirror dark matter ( 12 ) and hidden sector dark 
matter ( 12 - 16 ) all predict anisotropic scattering with ~ Ibarn/GeV = 0.6cm^/g, 

similar to nuclear cross-sections in the standard model. Note that couplings within the dark 
sector can be many orders of magnitude larger than those between dark matter and standard 
model particles, which is at most of order picobams ( 17 ). 

In terrestrial collider experiments, the forces acting on particles can be inferred from the 
trajectory and quantity of emerging material. Collisions between galaxy clusters, which contain 
dark matter, provide similar tests for dark sector forces. If dark matter’s particle interactions 
are frequent but exchange little momentum (via a light mediator particle that produces a long- 
ranged force and anisotropic scattering), the dark matter will be decelerated by an additional 
drag force. If the interactions are rare but exchange a lot of momentum (via a massive mediator 
that produces a short-ranged force and isotropic scattering), dark matter will tend to be scattered 
away and lost ( 11 , 18 , 19 ). 

The dynamics of colliding dark matter can be calibrated against that of accompanying stan¬ 
dard model particles. The stars that reside within galaxies, which are visible in a smoothed map 
of their optical emission, have effectively zero cross-section because they are separated by such 
vast distances that they very rarely collide. The diffuse gas between galaxies, which is visible 
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in X-ray emission, has a large electroweak cross-section; it is decelerated and most is eventu¬ 
ally stripped away by ram pressure (20). Dark matter, which can be located via gravitational 
lensing (21), behaves somewhere on this continuum (Fig. 1). 

The tightest observational constraints on dark matter’s interaction cross-section come from 
its behavior in the giant ‘bullet cluster’ collision 1E0657-558 (22). A test for drag yields 
o'DM/'m- < 1.25cm^/g (68% CL), and a test for mass loss yields < 0.7cm^/g (68% 

CL) (18). Half a dozen more galaxy cluster collisions have since been discovered, but no tighter 
constraints have been drawn. This is because the analysis of any individual system is fundamen¬ 
tally limited by uncertainty in the 3D collision geometry (the angle of the motion with respect 
to our line of sight, the impact parameter, and the impact velocity) or the original mass of the 
clusters. 

The same dynamical effects are also predicted by simulations in collisions between low- 
mass systems (11). Observations of low-mass systems produce noisier estimates of their mass 
and position (23-25), but galaxy clusters continually grow through ubiquitous ‘minor mergers’, 
and statistical uncertainty can be decreased by building a potentially very large sample (26,27). 
Lurthermore, we have developed a statistical model to measure dark matter drag from many 
noisy observations, within which the relative trajectories of galaxies, gas, and dark matter can 
be combined in a way that eliminates dependence upon 3D orientation and the time since the 
collision (28). 

We have studied all galaxy clusters for which optical imaging exists in the Hubble Space 
Telescope (Advanced Camera for Surveys) data archive (29) and X-ray imaging exists in the 
Chandra Observatory data archive (30). We select only those clusters containing more than 
one component of spatially extended X-ray emission. Our search yields 30 systems, mostly 
between redshift 0.2 < z < 0.6 plus two at z > 0.8, containing 72 pieces of substructure in 
total (Table SI). In every piece of substructure, we measure the distance from the galaxies to 
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the gas (5 sg- Assuming this lag defines the direction of motion, we then measure the parallel (5si 
and perpendicular (5 di distance from the galaxies to the lensing mass (Fig. 2). 

We first test the null hypothesis that there is no dark matter in our sample of clusters (a 
similar experiment was first carried out on the Bullet Cluster, finding a 3.4 and Sa detection 
(31)). Observations that do not presuppose the existence of dark matter (32) show that IO^^Mq 
clusters contain only 3.2% of their mass in the form of stars. We compensate for this mass, 
which pulls the lensing signal towards the stars and raised (5 gi by an amount typically 0.78 ± 
0.30 kpc (computed using the known distances to the stars ^sg; see Materials and Methods). 
The null hypothesis is that the remaining mass must be in the gas. However, we observe a 
spatial offset between that is far from the expected overlap, even in the presence of combined 
noise from our gravitational lensing and X-ray observations (Fig. 3A). A Kolmogorov-Smirnov 
test indicates that the observed offsets between gas and mass are inconsistent with the null 
hypothesis at 1.6a, a p-value of 3 x (without compensation for the mass of stars, this 

is 1 .la). This test thus provides direct evidence for a dominant component of matter in the 
clusters that is not accounted for by the luminous components. 

Having reaffirmed the existence of dark matter, we attempt to measure any additional drag 
force acting upon it, caused by long-range self-interactions. We measure the spatial offset of 
dark matter behind the stars, compensating as before for the 16% of mass in the gas (33) by 
subtracting a small amount from (5si (on average 4.3 ± 1.6 kpc). We measure a mean dark 
matter lag of ((5si) = —5.8 ± 8.2 kpc in the direction of motion (Fig. 3B), and ((5 di) = 1-8 ± 
7.0 kpc perpendicularly. The latter is useful as a control test: symmetry demands that it must 
be consistent with zero in the absence of systematics. We also use its scatter as one estimate of 
observational error in the other offsets. 

We interpret the lag through a model (28) of dark matter’s optical depth (similarly to pre¬ 
vious studies (19, 23)). Gravitational forces act to keep gas, dark matter and galaxies aligned. 
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while any extra drag force on dark matter induce a fractional lag 



where Ugai is the interaction cross-section of the galaxies, coefficient B encodes the relative 
behavior of dark matter and gas, and a*/m is the characteristic cross-section at which a halo of 
given geometry becomes optically thick. We assume that stars do not interact, so agai 0. To 
ensure conservative limits on a-Qu/rn, we also assume 5^1 and marginalize over a*/m ^ 
6.5 ± 3cm^/g, propagating this broad uncertainty to our final constraints (see Materials and 
Methods). Adopting the dimensionless ratio /3 brings two advantages. First, it removes de¬ 
pendence on the angle of the collision with respect to the line of sight. Second, it represents 
a physical quantity that is expected to be the same for every merger configuration, so mea¬ 
surements from the different systems can be simply averaged (with appropriate noise weight¬ 
ing, although in practice, the constraining power from weak lensing-only measurements comes 
roughly equally from all the systems). 

Combining measurements of all the colliding systems, we measure a fractional lag of dark 
matter relative to gas {f3) = —0.04±0.07 (68% CL). Interpreting this through our model implies 
that dark matter’s momentum transfer cross-section is anu/'m = —0.25lo;43 cm^/g (68% CL, 
two-tailed), or aBM/m < 0.47cm^/g (95%CL, one-tailed); the full PDF is shown in Fig. 4. 
This result rules out parts of model space of hidden sector dark matter models e.g. (12,13,15,16) 
that predict ~ 0.6 cm^/g on cluster scales through a long-range force. The control test 

found {/3±) = {Sdi/6sg) = —0.06 ± 0.07 (68% CL), consistent with zero as expected. This 
inherently statistical technique can be readily expanded to incorporate much larger samples 
from future all-sky surveys. Equivalent measurements of mass loss during collisions could also 
test dark sector models with isotropic scattering. Combining observations, these astrophysically 
large particle colliders have potential to measure dark matter’s full differential scattering cross- 
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Hot, diffuse gas 

visible in X-rays 


Dark matter 



Figure 1: Cartoon showing the three eomponents in eaeh pieee of substrueture, and their relative 
offsets, illustrated by blaek lines. The three eomponents remain within a eommon gravitational 
potential, but their eentroids beeome offset due to the different forees aeting on them, plus 
measurement noise. We assume the direetion of motion to be defined by the veetor from the 
diffuse, mainly hydrogen gas (which is stripped by ram pressure) to the galaxies (for which 
interaction is a rare event). We then measure the lag from the galaxies to the gas and to the 
dark matter in a parallel ^si and perpendicular ^di direction. 
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Figure 2: Observed eonfigurations of the three eomponents in the 30 systems studied. The 
baekground shows the HST image, with eontours showing the distribution of galaxies (green), 
gas (red) and total mass, whieh is dominated by dark matter (blue). 
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Observed offset between various components of substructure [kpc] 

Figure 3: Observed offsets between the three eomponents of 72 pieees of substrueture. Offsets 
(5si and ^gi inelude eorreetions aeeounting for the fact that gravitational lensing measures the 
total mass, not just that of dark matter. (A) The observed offset between gas and mass, in the 
direction of motion. The smooth curve shows the distribution expected if dark matter does not 
exist; this hypothesis is inconsistent with the data at 1.6a statistical significance. (B) Observed 
offsets from galaxies to other components. The fractional offset of dark matter towards the gas, 
<5 si/<5sg, is used to measure the drag force acting on the dark matter. 
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Dark matter self-interaction cross section, acy^lm [cm^/g] 


Figure 4: Constraints on the self-interaetion eross-seetion of dark matter. These are derived 
from the separations /3 = 5 si/<^sg> assuming a dynamieal model to eompare the forees aeting 
on dark matter and standard model partieles {28). The hatehed region denotes 68% eonfidenee 
limits, to be eompared to the 68% eonfidenee upper limits from previous studies of the most 
eonstraining individual elusters in blue. Note that the tightest previous eonstraint is derived 
from a measurement of dark matter mass loss, whieh is sensitive to short range self-interaetion 
forees; all other eonstraints are measurements of a drag foree aeting on dark matter, eaused by 
long range self-interaetions. 
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Materials and methods 

We followed overall proeedures that we developed in blind tests on moek data (24), usually exploiting 
algorithms for high preeision measurement that had been developed, ealibrated and verified elsewhere. 
However, several eustom adaptations were required to analyze the heterogeneous data from the Hubble 
Spaee Teleseope (HST) and Chandra X-ray Observatory arehives (Table S1 lists all the observed systems 
(34-37), and Figure S2 shows the offsets measured in eaeh). 

Here we deseribe the methods we used to eombine observations with different exposure times, fil¬ 
ters, epoehs and orientations - starting from the raw data and performing a full reduetion to maximise 
data quality. To eonvert angular distanees into physieal distanees, we assume a eosmologieal model 
derived from measurements of the Cosmic Microwave Background (33), Om = 0.31, Oa = 0.69, 
Ho = 67km/s/Mpc. 


Position of gas, seen in X-ray emission 

We downloaded the raw event 1 files for all observations. To process fhese dafa, we used CIAO fools 
version 4.5, sfarfing wifh basic reduction and calibrafion using fhe CIAO repro fool. In our analysis, if is 
parficularly imporfanf fo remove emission from poinf sources, and prevenf X-ray brighf Acfive Galacfic 
Nuclei af fhe cenfers of clusters from biasing our posifion measuremenfs. We Iherefore made a firsf pass 
af removing poinf sources using celldetect. We fhen filtered each even! fable for any pofenfial spurious 
evenfs such as solar flares by clipping fhe fable af fhe Aa below fhe mean flux level. 

Having cleaned each exposure, we combined fhem using fhe merge-obs scripl from CIAO fools info 
a single exposure map correcfed flux image, producing along wifh if exposure maps for each observafion 
and fhe sfacked image. We modeled fhe Chandra PSF af each posifion fhroughouf fhe field, we creafed 
individual maps using mkpsfmap for each exposure af an effecfive energy of 1 keV, fhen combined each 
model weighfing fhem by fheir respecfive exposure map. Figure S3 shows an example of fhe PSF map 
used for cluster A520. 

To make a second pass fo idenlify poinf sources, we passed fhe sfacked image and PSF model fhrough 
CIAO wavdetect, a wavelef smoofhing algorifhm fhaf employs a ‘Mexican haf’ filter on a range of scales. 
This esfimafes fhe frue size of each source, correcting for fhe size of fhe PSF. We used fhe smallesf 
scales for fhe wavelef radii (1, 2 pixels) fo idenlify poinf sources, and combined fhe larger scales (4, 
8, 16, and 32 pixels) info a denoised version of fhe final image. We finally inspecfed every image by 
eye for any remaining poinf sources. We found lhal Ibis double filler melhod proved very successful af 
removing poinf sources, wifh only AGN af fhe edge of fhe clusler remaining unfiagged. Allhough fheir 
emission has extended wings, fhe clusler is usually in fhe center of fhe poinling, resulling in minimal 
contamination. 

Finally, we measured the position of coherent substructure in the X-ray emission using SExtrac- 
tor (38). This calculates positions from the first order moments of the light profile, which means lhal fhe 
refurned posifion does nol always coincide exaclly wifh fhe brighlesl pixel. SExtractor does nol reporl 
reliable errors in fhe positions bul, since fhe dominanl conlribulion of varialion is fhe size of fhe smoolh- 
ing kernel, we can estimate fhe robuslness of our measuremenfs by smoofhing fhe image using differenl 
scales in wavdetect, and measure fhe rms across differenl scales. On average we found fhe rms error fo 
be 4 arcseconds (roughly 30kpc af redshifl ^=0.4). 
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Cluster 

RA (deg) 

DEC (deg) 

z 

ACS Filter 

ACS (s) 

Chandra (ks) 

1E0657 

104.612 

-55.9477 

0.296 

F814W_F775W 

15094.0 

597.39 

A1758 

203.194 

50.5426 

0.2792 

F814W 

10000.0 

216.00 

A209 

22.9728 

-13.6127 

0.206 

F814W 

4040.00 

24.09 

A2146 

239.007 

66.3725 

0.234 

F814W 

9233.00 

84.08 

A2163 

243.938 

-6.14690 

0.203 

F814W 

9192.00 

444.59 

A2744 

3.58210 

-30.3898 

0.308 

F814W 

11980.0 

133.12 

A370 

39.9627 

-1.58000 

0.373 

F814W 

3840.00 

86.81 

A520 

73.5395 

2.93110 

0.202 

F814W 

18320.0 

426.01 

A781 

140.149 

30.4927 

0.298 

F814W 

1620.00 

49.54 

ACTCLJ0102 

15.7277 

-49.2560 

0.87 

F814W 

1916.00 

359.16 

DLSCLJ0916 

139.046 

29.8450 

0.5343 

F814W 

9894.00 

41.28 

MACSJ0025 

6.37460 

-12.3818 

0.5843 

F814W 

4200.00 

168.61 

MACSJ0152 

28.1473 

-28.8944 

0.341 

F606W 

1200.00 

20.04 

MACSJ0358 

59.7174 

-29.9320 

0.428 

F814W 

4620.00 

65.74 

MACSJ0416 

64.0392 

-24.0735 

0.42 

F814W 

4037.00 

57.50 

MACSJ0417 

64.3926 

-11.9111 

0.443 

F814W 

1910.00 

95.92 

MACSJ0553 

88.3494 

-33.7117 

0.407 

F814W 

4572.00 

88.74 

MACSJ0717 

109.389 

37.7528 

0.5458 

F814W 

8893.00 

83.22 

MACSJ1006 

151.730 

32.0198 

0.359 

F814W 

1440.00 

13.30 

MACSJ1226 

186.694 

21.8673 

0.37 

F814W 

5520.00 

153.81 

MACSJ1354 

208.635 

77.2528 

0.3967 

F814W 

1200.00 

35.46 

MACSJ1731 

262.913 

22.8660 

0.389 

F814W 

1440.00 

22.28 

MACSJ2243 

340.837 

-9.58910 

0.447 

F606W 

1200.00 

21.88 

MS1054 

164.245 

-3.62000 

0.826 

F606W 

8100.00 

89.51 

RXCJ0105 

16.4096 

-24.6801 

0.23 

F606W 

1200.00 

21.97 

RXCJ0638 

99.6953 

-53.9735 

0.1658 

F606W 

1200.00 

21.78 

RXJIOOO 

150.132 

44.1491 

0.154 

F606W 

1200.00 

20.66 

SPTCL2332 

352.959 

-50.8642 

0.5707 

F606W 

7680.00 

39.9 

ZWCL1234 

189.045 

28.9929 

0.2214 

F814W 

27632.0 

51.75 

ZWCL1358 

209.951 

62.5163 

0.329 

F850LP 

13692.0 

63.10 
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Figure SI: The full sample of 30 merging complexes, and their locations on the sky. The 
columns show, from left to right: the name of the cluster, its right ascension, declination, and 
redshift, the HST/ACS filter used for our lensing analysis, and the total exposure time for that 
particular filter, and the (cleaned) exposure time of the Chandra X-ray image. 
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Figure S2: Observed offsets between galaxies, gas and dark matter in 72 eomponents of sub¬ 
structure. In each case, the green triangle, at the centre of the coordinate system, denotes the 
position of the galaxies. The separation between galaxies and gas, is shown in red. The 
separation of the dark matter with respect to the galaxies, projected onto the SG vector, (isi, is 
shown in blue. The error bars show the locally estimated la errors. 
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Figure S3: An example model of the size of the Chandra X-ray teleseope’s Point Spread Fune- 
tion (PSF). The model PSF is used to identify and remove point sourees, e.g. Aetive Galaetie 
Nuelei - and to thereby identify extended X-ray emission from hot gas within the eluster. The 
image shows a combined, exposure map weighted, PSF map stacked for the various observa¬ 
tions of galaxy cluster A520. 
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Position of galaxies, seen in optical emission 

We searched the HST archive for data acquired with the Advanced Camera for Surveys (ACS) instrument, 
which has the largest field of view. We considered only filters F606W, F814W and F850LR, whose high 
throughput ensures deep imagining, and whose red wavelengths ensure (a) that the optical emission 
samples the old stars that dominate the mass content of these systems and (b) a high density of high 
redshift galaxies visible behind the cluster, to provide sufficient lensing signal. Some clusters had been 
observed in more than one wavelength band. We used only a single band for all the clusters to further 
homogenize the data, but have compared a subset of our results in different bands to check for systematic 
errors. For our main analysis, we selected the broad F814W band, unless there are significantly more 
exposures in another. 

We corrected the raw, pixellated data for charge transfer inefficiency (39), fhen performed basic data 
reduction and calibration using the standard Calacs pipeline. We used tweakReg to orient and align 
individual exposures, then stacked them using MultiDrizzle (40) with a Gaussian convolution kernel and 
PIXFRAC=0.8 (41) to produce a deep, mosaicked image with a pixel scale of 0.03 arcseconds. In the 
process, MultiDrizzle also output a reoriented image of each individual exposure, which we used for 
star/galaxy identification and PSF estimation. 

We estimated the distribution of mass in galaxies via the proxy of the light emitted by their stars. In 
our single-band imaging, we were able to identify and mask foreground stars in the Milky Way (which 
appear pointlike), but assumed any foreground or background galaxies to be randomly positioned and 
thus merely add shot noise to our measurements. We smoothed the masked image using wavdetect, and 
measured the position of coherent substructure using SExtractor (38). This calculates positions from the 
first order moments of the light profile, which means fhat the returned position does not always coincide 
exactly with the brightest pixel. SExtractor does not report reliable errors in the positions. However, 
since the dominant contribution of noise is inclusion or omission of galaxies inside the smoothing ker¬ 
nel, we estimated the robustness of our measurements by smoothing the image using different scales in 
wavdetect, and compared the resulting positions. On average, we found an rms error in the position of 
the extracted halos of 0.6 arcseconds (roughly 4.5 kpc at redshift z=0.4). 

We also tried two other ways to quantify the position of the galaxies. First, we measured the 
smoothed distribution of galaxies in the image, with all galaxies weighted equally (this represents the 
opposite - and least realistic - assumption of galaxies’ mass/light ratio). To do this in practice, we 
passed the galaxy catalogue through the X-ray data reduction pipeline, as if each galaxy were a single 
X-ray photon. This created a smoothed image, in which we identified substructure using SExtractor. 
Since the same galaxies contributed both to the flux-weighted and galaxy-weighted positions, the two 
measurements are correlated. We measure the uncertainty on the galaxy weighted positions to be 5 kpc, 
about the same as the flux-weighted positions. We obtain consistent values of (/)) = 0.054 ± 0.062 (68% 
CL) and conclude that (Tdm/w^ = 0-36lo45 cm^/g (68% CL, two-tailed). Second, we tried identifying 
the position of the ‘Brightest Group Galaxy’ (BCG), since its formal error is small, and it has proved 
optimal in studies of isolated groups (42). In merging systems however, the brightest nearby galaxy is 
frequently unassociated with the infalling group (43). Accounting for our observed 1.7 ± 0.9 arcsecond 
offset to any brighter galaxy within 25 arcseconds of X-ray emission (the search region that will be used 
to identify gravitational lensing signals), again yields a consistent constraint on a-QM/'m,, but with much 
larger final error. 
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Position of dark matter, measured via weak gravitational lensing 

We measured the elliptieities of galaxies in HST images using the RRG method (44). This eorreets 
galaxies’ Gaussian-weighted moments for eonvolution with the Point Spread Function (PSF), to measure 
the shear 71 ( 72 ) corresponding to elongations along (at 45 degrees to) the x axis. This method has been 
empirically calibrated on simulated HST imaging in which the true shear is known (45), applying a 
multiplicative correction of (m) = —3.0 x 10“^ and a additive bias of (c) = —2.1 x 10“^. 

HST’s PSF varies across the field of view and, because thermal variations change the telescope’s 
focus, at different epochs. Modelling the net PSF in our stacked images therefore required a flexible 
procedure. We first identified stars in the deep, stacked image using their locus in size-magnitude space. 
We then measured the ellipticity of each star in individual exposures. By comparing these to TinyTim {46) 
models of the HST PSF (created by raytracing through the telescope at different focus positions but at 
the appropriate wavelengths for the band), we determined the focus position for each exposure. We then 
interpolated (second and fourth shape moments of) the TinyTim PSF model to the position of the galaxies, 
rotating into the reference frame of the MultiDrizzle mosaic. We then summed the PSF moments from 
each exposure in which a galaxy was observed. Figure S4 shows an example of the final PSF model for 
one cluster. 

We measured the shear of all galaxies that appear in 3 or more exposures, with a combined signal-to- 
noise in the stacked image > 4.4 and size >0.1 arcseconds. These cuts (45) remove noisy measurements 
at the edges of the field or in the gaps between detectors. We also masked out galaxies that lie near bright 
stars or large galaxies, whose shapes appear biased. Figure S5 compares shear catalogues for a single 
cluster, derived from independent analyses of data in the F814W and F606W bands. There is the expected 
level of scatter between the two measurements - but, most importantly, there is no detectable bias. 

We reconstructed the distribution of mass in the clusters using the parametric model-fitting algo¬ 
rithm Lenstool {47). Using Bayesian likelihood minimization, Lenstool simultaneously fits multiple mass 
haloes to an observed shear field, with the position and shape of each halo described by the NFW {48) 
density profile. This is an efficient technique to record a unique position for each halo, marginalizing 
over nuisance parameters that include mass and morphology, that are not of direct interest to our study. 
Assuming this density profile does not bias measurements of the position of halos within current statis¬ 
tical limits (24). Lenstool requires positional priors to be defined in which it searches for the lensing 
signal. Except in a few well-studied systems (where we use the extra information), we obtained an initial 
lensing model using one prior search radius centered on each gas position and large enough to incorpo¬ 
rate any nearby groups of galaxies. Following this scheme, we used an automated procedure to identify 
and associate the mutually closest galaxy, gas and lensing signals into systems of three mass components. 
In all systems, we then modeled the lensing system a final time, adopting priors centred on the galaxy 
position (we redid this step when trying different position estimators for the galaxies). Henceforth, we 
could center the coordinate system for each combined system of galaxies, gas and dark matter on the 
galaxies, to avoid prior bias in the Bayesian fits. 

Lenstool samples the posterior surface in two ways. To obtain the best fitting position, we iterated 
to the best-fit solution with a converging MCMC step size, using ten simultaneous sampling chains to 
avoid local maxima. To sample the entire posterior surface (whose width quantifies uncertainty on model 
parameters), we then reran the algorithm with a fixed step size. The la error on position was on average 
11.4 arcseconds (roughly 60 kpc at redshift z=0.4). As a sanity check we compare our measured centroids 
to those systems included in previous studies. Our statistical uncertainty is sometimes larger because we 
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Figure S4: An example model of the Point Spread Funetion (PSF) of the Hubble Spaee Tele- 
seope/Advaneed Camera for Surveys (HST/ACS). Eaeh tiek mark represents the elliptieity of 
the PSF at that partieular position in the HST field. Its orientation shows the PSF’s major axis 
and its length shows the elliptieity; a dot would indieate a eireular PSF The PSF tends to be 
highly elliptieal near the edge of the field and more eireular in the eentre. Tiek marks are plot¬ 
ted at the position of every “deteeted” souree. The mosaie pattern of dithered exposures ean be 
seen: noisier regions with fewer exposures eontain more spurious sourees, whieh are removed 
during analysis (but are shown here for elarity). The example shown is for observations of 
galaxy eluster MACSJ0416. 
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Figure S5: A comparison of the gravitational lensing shears measured independently behind a 
single cluster, in two different HST filters. The top (bottom) panel shows the difference between 
7i ( 72 ) for each galaxy, which traces to elongations along (at 45 degrees to) the x axis. We find 
scatter as expected due to observational noise, but no systematic bias. 
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use only weak gravitational lensing, but we find no evidenee for any bias. For example, our measured 
positions in the ‘bullet eluster’ lie within one standard deviation of those reported in {31). 


Positional offsets between components 

When assigning different mass eomponents to one another, for almost all the elusters, we used an auto¬ 
mated matehing algorithm to assoeiate the nearest elumps of dark matter, gas and stars. This was made 
robust by performing the matehing in both direetions (e.g., dark matter to stars, and stars to dark matter). 
In a few eases where detailed analyses of individual systems were available in the literature (for exam¬ 
ple, using strong lensing. X-ray shoeks, optieal speetroscopy or imaging additional bands, which were 
outside the scope of our work), we inserted that prior information by hand during association. This was 
most useful in systems A520 and A2744. As a further test, we carry out a jackknife test to ensure that the 
association does not effect the overall constraints, and moreover, no single cluster dominates the result. 
We find no evidence for such an effecf, and derive consisfenf error bars of AcjDM,jK/?7r = ±0.5cm^/g, 
furfher supporfing fhe error bars quofed in our final resulf. 

We drew an offsel vector (isG in angle befween fhe observed position of fhe gas and galaxies, which 
we took fo define fhe sysfem’s direcfion of mofion. We fhen measured fhe position of fhe fofal mass along 
fhaf vector and (in a righf handed coordinate system) perpendicular fo if, defining offsef vecfors Jsi, (Jgg 
and (5 di from fhe infersecfion poinf I of fhese vecfors. 

Gravifafional lensing measures fhe posifion of lofal mass, rafher fhan fhaf of Jusf dark mailer. We 
corrected fhe measured offsels Jsi and (5 gi for fhe conlribulion from fhe nexl mosl massive componenl. 
To calibrale Ihis correcfion, we analysed mock lensing dala from a dominanl mass componenl (wilh 
an NFW {48) profile) plus a less massive componenl al some offsel 5. The correclions were always 
small bul, for a subdominanf componenl wilh fhe same profile, normalised fo conlain a fraclion / of fhe 
lofal mass, we found fhaf fhe lensing position is pulled by an amounf , and we corrected 

for fhaf. If we do nol calibrale for fhe exlra pull of gas on fhe lensing peak we infer an upper limil of 
< 0.54cm^/g ( 68 % CL, one-lailed). 

To lesl fhe hypolhesis fhaf dark mailer does no! exisl, we required a model of fhe (Jgi dala expecfed 
if Ihis were Irue. To generale fhaf model, we assumed fhaf fhe Irue posilions of fhe X-ray and lensing 
signals coincided, bul fhaf fhe observed posilions were offsel by a random amounl determined by Ihe 
appropriate level of noise in each (see above). We calculated Ihe 2D offsel, Ihen projected Ihis onto Ihe 
direction to Ihe slars, which is also selected al random. We could have slighlly increased Ihe model (5 gi 
offsel to accounl for Ihe mass in slars (Ihe increase musl be positive because Ihe vector (5 sg is defined 
from Ihe galaxies to Ihe gas). However, if is belter to instead decrease Ihe observed (5 gi offsel. The Iwo 
approaches are equivalenl in principle, bul Ihe laller allowed information to be added to our analysis 
because Ihe absolute value of (5sG was known in each system. When comparing Ihe model and observed 
(5 gi offsels via a Kolmogorov-Smirnov lesl (in which we computed critical values using a Taylor series), 
we also used Ihe errors on itgi determined for each system individually. 

When measuring Ihe interaction cross-section of dark mailer, we converted offsel measuremenls in 
arcseconds to physical unils of kpc (using a standard cosmological model, which assumes dark mailer 
exisls). This enabled a more detailed comparison of Ihe offsels belween differenl systems. The (nois¬ 
ily determined) error estimates of offsels in a few systems were anomalously low, and likely smaller 
lhan Ihe uncertainly in our knowledge of Ihe merging configuration. To more robuslly quantify Ihe to¬ 
tal uncertainly of offsels (which should include observational noise plus Ihe possibility of componenl 
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misidentification and merging irregularities), we empirieally exploited the eontrol test 5 di, whieh has 
an rms variation between systems ctdi = 60kpe. This value is eonsistent with most of the individually 
measured errors, but more robust. We therefore adopted it globally as the error on every measurement 
of (5 di and (5si, resealing to a value in areseeonds at the redshift of eaeh system. Errors in 5 sg must be 
smaller than this, beeause they do not involve observational noise in the lensing position. However, they 
also inelude the possibility of eomponent misidentiheation, whieh is best estimated through this global 
approaeh. We therefore adopted the eonservative approaeh of also assigning this value as the error on 
every measurement of (5 sg- Thus we set csg = o'si = cjdi = 60kpe. To eombine our measurements of 
/5 = 5 si/<^SG and I3± = (5di/(^SG from individual systems, we multiplied their posterior probabilities (ap¬ 
proximated as a normal distribution even though it is a Cauehy distribution, but with a width determined 
by propagating errors on the individual offsets). 


Interpreting positional offsets as an interaction cross-section 

Similarly to previous studies of the eross-seetion of dark matter {19,23), we interpreted observations of 
offset dark matter in terms of its optieal depth for interaetions. However, we have developed a more 
sophistieated model {28) intended to take into aeeount the 3D and time-varying trajeetories of infalling 
halos. First, ealeulating the dimensionless ratio j3 = 5si/()sG removes dependenee on the angle of the 
eollision with respeet to the line of sight. Furthermore, a set of analytie assumptions suggests that j3 
is a physieally meaningful quantity that should be the same for every system. The main assumption of 
quasi-steady state equilibrium is reasonable for the deteetable systems in our sample, but eaution would 
be needed to interpret dark matter substrueture that had passed direetly through the eluster eore (and 
had its gas stripped) or substrueture on a radial orbit eaught at the brief moment of turnaround (this is 
a negligible fraetion in our moek data). The model also ineorporates the results of simulations {11) in 
whieh dark seetor interaetions that are frequent but exehange little momentum (e.g. via a light mediator 
partiele that produees a long-ranged foree and anisotropie seattering) produee a drag foree and separate 
dark matter from the stars. On the other hand, simulations of ‘billiard ball’ interaetions that are rare 
but exehange a lot of momentum (e.g. via a massive mediator that produees a short-ranged foree and 
isotropie seattering) tend to seatter dark matter away from a system and produee mass loss {11,18,19). 
However, we note that the ref. {18) also reports an unexpeeted small separation between galaxies and 
dark matter after billiard ball seattering. In this paper, we explieitly follow the preseription in {11). 

Aeeording to our model of dark matter dynamies (see equation 33 of ref. {28)), the offset of dark 
matter from galaxies, ealibrated against the offset of gas, is 

= (SI) 

Sinee the gaps between galaxies are vast eompared to their size, they interaet very rarely, so we assumed 
that (Tgai ~ 0. If this assumption were wrong, or in the presenee of observational noise, our analysis 
ean therefore produee negative values of udm/w.. Our quoted errors inelude observational errors, pair- 
assignment errors, and model parameter errors. 

The value of a*jm depends upon the geometrieal properties of the dark matter halo, but is pro¬ 
portional to its mass and inversely proportional to its eross-seetional area. For our set of merging 
systems, we eonservatively adopted a a*jm = 6.5 ± 3cm^/g by assuming the system masses are 


11 



log(M 2 oo/-^ 0 ) = 14 ± 1, with NFW density profiles and eoneentration varying with mass as ob¬ 
served in numerieal simulations (49). By assuming a eonservative range in halo masses we propagate 
a mueh larger error in a* jm than one would expeet if we were to measure the true values. We then 
analytieally marginalized over a* jm, propagating the uneertainty through to our final eonsfrainfs. The 
lop panel of Figure S 6 shows fhe values of a* jm instead assuming differenl, fixed system masses; Ihe 
bottom panel shows fhe effeel on udm/?^- The inferred estimate of aDM/^n is broadly insensitive lo 
a*jm, varying from aYyyijvn = —0.23 ± 0.60cm^/g for an assumed halo of M 200 = lo 

= —0.1 ± 0.28 cm^/g for a halo of M 200 = IO^^Mq. 

The relalive behaviour of gas and dark matter was eompared Ihrough a ratio in fhe prefaelor 

^ _ C'pM^DMll^gasPDM 

C*gas ^gas Pgas 

where C, A and M are fhe drag eoeffieienl, size and mass of fhe merging halo, and p is fhe density 
of malerial Ihrough whieh if is moving. We assumed a eonservalive lower limil of i? 1, leading lo a 
eonservative upper limil on our eonslrainls on (Tdm/?71- 

The lirsl requiremenl lo have ensured a eonservative Irealmenl is lhal Ihe infalling subslruelure’s gas 
envelope is smaller lhan ils dark matter envelope, ylgas < ^dm- This is generieally Irue of isolated 
slruelures in numerieal simulations and, as gas is slripped during Ihe eollision, il will beeome smaller 
still. The geomelrie size of a gas halo also depends upon ils lemperalure - and hoi gas may be more easily 
slripped lhan eold gas. To lesl whelher lhal has an slatislieally signilieanl effeel, we measured Ihe X-ray 
lemperalure of eaeh observed infalling system, and separately analyzed Ihe hotter and eooler halves of 
our sample. As shown in Figure S7, Ihe resulls for eaeh half remain eonsislenl, wilh error bars larger 
by approximately \/2. For Ihe hotter sample (T > 8 keV), we found OT^yilm, = —0.10 ± 0.58cm^/g 
and for Ihe eooler sample (T < 8 keV) we found (Tdm/™. = —0.50 ± 0.64 cm^/g. Allhough Ihere was 
marginal evidenee lhal hoi gas is more easily slripped lhan eold gas, whieh eould be investigated wilh a 
mueh larger sample, our eonelusions remain unaffeeled wilhin eurrenl slalislieal preeision. 

The seeond requiremenl lo have ensured a eonservative Irealmenl is lhal Ihe gas fraelion in Ihe 
medium Ihrough whieh Ihe bullel is Iraveling, /ga^ = 9gas/ Pdm is less lhan lhal of Ihe infalling slrue- 
lure Mgas/MDM- We assumed lhal, overall, infalling slruelure eonlains Ihe universal fraelion = 

0.184 (33), and we measured /gas in moek dala realised from eosmologieal simulations of slruelure 
formation (50). The mean /gas over all simulations (Ihe solid line in Figure S 8 ) is lower lhan Ihe uni¬ 
versal fraelion, and is indeed eonslanl (wilhin 10 %) al Ihe radii of observable subslruelures (poinls in 
Figure S 8 ). These eonelusions from simulations are eonsislenl wilh deep X-ray observations of galaxy 
eluslers, e.g. (5i). 
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Figure S6: The sensitivity of measurements of dark matter’s self-interaetion eross-section to the 
model parameter a*/m. This parameter is the eharaeteristic value of eross-seetion at whieh an 
appropriately-sized eloud of standard model partieles beeomes optieally thiek. The top panel 
shows the value of a*/m for different various substructure masses, assuming an NFW mass pro¬ 
file and a mass-concentration relation from cosmological simulations (49). The bottom panel 
demonstrates how a few of those values affect our measurement of the cross-section. The re¬ 
sulting variation is sub-dominant to statistical error in our sample of clusters. We adopted a 
value of a*/m = 6.5 ± 3cm^/g, corresponding to dark matter halos of M = and 

propagated the uncertainty through to our final constraints. 
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Figure S7: The sensitivity of measurements of dark matter’s self-interaetion eross-seetion to 
the temperature of the gas against whieh dark matter’s trajeetory is ealibrated. We measured 
the gas temperature from the X-ray speetra of our 72 systems, and split the sample in two: blue 
data show substructures with gas temperature < 8 keV, and red data show substructures with 
gas temperature > 8 keV. The constraining power of each sample is approximately \/2 less than 
that of the full sample, shown in grey, and no statistically significant difference is measured. 
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Figure S8: The sensitivity of measurements of dark matter’s self-interaction cross-section to the 
density of gas through which it is moving. The plot shows the gas fraction /gas = Pgas/poM 
in simulated galaxy clusters (50), as a function of clustercentric radius. The solid line shows 
the average /gas over 16 clusters, with the la error on the mean given in grey. Triangles show 
the measured /gas at the radius of substructures observable in mock 2D realisations of the 3D 
simulations (only the inner ~ 60% lie inside the HST field of view at the redshifts of the 
observed systems). Our interpretation of the dark matter and gas trajectories as an interaction 
cross-section, assumes that these are lower than the universal fraction Ob/Od = 0.184 (33). 
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